Exploiting soliton decay and phase fluctuations in atom chip interferometry of 

Bose-Einstein condensates 
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We show that the decay of a soliton into vortices provides a mechanism for measuring the initial 
phase difference between two merging Bose-Einstein condensates. At very low temperatures, the 
mechanism is resonant, operating only when the clouds start in anti-phase. But at higher tem- 
peratures, phase fluctuations trigger vortex production over a wide range of initial relative phase, 
as observed in recent experiments at MIT. Choosing the merge time to maximize the number of 
vortices created makes the interferometer highly sensitive to spatially varying phase patterns and 
hence atomic movement. 

PACS numbers: 03.75.Kk, 03.75.Lm 
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Recent experiments [l[ have successfully performed in- 
terferometry of two Bose-Einstein condensates (BECs) in 
an atom chip trap. The two BECs were merged over a 
precisely controlled timcscalc and then allowed to expand 
freely. The size of the cloud after expansion was found to 
increase with the initial phase difference, A, between the 
two BECs, being maximal when A = n. Consequently, 
measurements of the cloud size following expansion can 
be used to determine A, and hence form the basis of an 
atom interferometer. However, the physical mechanism 
that makes the expanded cloud size depend on A is un- 
known. As a result, many interesting questions remain 
unanswered, for example concerning the variation of ex- 
panded cloud size with merging time [l| . In addition, for 
the temperatures and extremely elongated geometry of 
the BECs studied in the experiments, phase fluctuations 
are predicted [3, La] , but it is unclear how this affects the 
operation of the interferometer. 

In this Letter, we show that vortex production is the 
physical mechanism that causes the size of the atom cloud 
after expansion to vary with A. When A = 7r, a soliton 
forms between the two merging BECs and then decays 
into vortices. The additional momentum associated with 
the vortices drives an enhanced spreading of the cloud 
when it is released from the chip trap. We investigate 
the role of the merging time on the creation of dynamical 
excitations and predict that the number of vortices pro- 
duced will be maximal for a merge time of ~ 5 ms. Sur- 
prisingly, the finite temperature phase fluctuations @, Q 
characteristic of elongated BECs can actually enhance 
the sensitivity of the interferometer. This is because, in 
the absence of phase fluctuations, vortex production is a 
sharply resonant process that only occurs when A = ir. 
However, the presence of phase fluctuations causes the 
number of vortices generated by the merging to change 
more smoothly with A. Consequently, the action of the 
interferometer can be controlled by varying the temper- 
ature and/or geometry of the BEC. 

In the experiments a 23 Na BEC containing N = 
4 x 10 5 atoms is held in an atom chip trap produced 



by three current-carrying wires and an external bias 
field. As expected for atom chips, the trap frequen- 
cies are very high perpendicular to the trapping wire 



2tt x 1000 rad s 1 ), but low along the wire 



(bj z = 2tt x 9 rad s _1 ), thus creating an exceedingly elon- 
gated cigar-shaped BEC of peak density n = 2.8 x 10 14 
cm~ 3 . For aspect ratios this extreme, finite temperature 
phase fluctuations prevent global coherence of the atom 
cloud, even at temperatures T << Tc, the transition 
temperature 1. In the limit of a truly one-dimensional 
gas, phase coherence is completely lost, even at T = |4|]. 
In the MIT experiments, T ~ 0.5 Tc, which is far higher 
than the threshold of 0.1 Tc above which phase fluctua- 
tions are expected [H. 

By applying an oscillating rf field from a fourth wire 
on the atom chip [B|, Q , the single- well potential is de- 
formed smoothly into a double-well potential, thus split- 
ting the BEC along its short (x) axis [Fig. Ufa)]. We 
assume that this process is performed adiabatically, and 
so begin our simulations with two equilibrium ground- 
state BECs, each containing 2 x 10 5 atoms, calculated by 
solving the three-dimensional Gross-Pitaevskii equation 
in imaginary time @, S| ■ I n the experiments [l[ , A is set 
by holding the atoms for a pre-determined time before be- 
ginning the merging process. Due to small asymmetries 
in the trapping potential, there is a reproducible differ- 
ence between the chemical potentials of the two clouds, 
which causes a non- vanishing and measurable phase evo- 
lution rate [(J. In our calculations, at the start of the 
merging process (time t = 0) we set the phase of the 
BEC to for x < and A for x > 0. The BECs are 
then merged by smoothly transforming the double-well 
potential back into a single-well potential over a known 
merging time r, as in experiment 

Throughout the merging process, we model the three- 
dimensional wavefunction by an expansion over plane- 
wave basis states Q, with a maximum cutoff wavevec- 
tor, k m , chosen to prevent Fourier aliasing 0]. In the ab- 
sence of fluctuations, this method is equivalent to solving 
the projected Gross-Pitaevskii equation [l(| • The Gross- 
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FIG. 2: Atom density profiles within the region of the y — 
plane (axes inset) shown by the dashed rectangle in Fig. 1(a), 
at t = t = 1 ms for A = (a) and tv (b), and at t = r = 50 
ms for A = (c) and it (d). Horizontal bar shows scale. 



FIG. 1: (a) Atom density profile of the two BECs in the y = 
plane (axes inset) at t — 0. (b)-(g) Density profiles within the 
region enclosed by the dashed rectangle in (a) at key stages of 
the merging process (r = 5 ms) calculated for A = at t — 3 
ms (b), 4 ms (c), 5 ms (d), and for A = n at t — 3 ms (e), 
4 ms (f), and 5 ms (g). Upper [lower] horizontal bars show 
scales in (a) [(b)-(g)]. 



Pitaevskii equation has been very successful in calculat- 
ing equilibrium states and modeling mean-field dynam- 
ics, but it fails to describe spontaneous processes, such 
as the formation of a scattering halo and the associated 
depletion of the BEC. We include such effects in our cal- 
culations by ad ding random complex noise to each plane 
wave mode [H, lllj The amplitude of these quantum 
fluctuations has a Gaussian distribution, with an aver- 
age value of half a particle. This appro ach is known as 
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the truncated Wigner method [§, ll|, 

In order to model the finite temperature phase fluc- 
tuations characteristic of elongated BECs, we add extra 
thermal fluctuations to low energy modes of the BEC 
below a small wavevector cutoff, k s = 1.3 x 10 6 m -1 , 
which is just sufficient to describe the overall curvature 
of the initial BEC wavefunction. To allow this extra 
noise to reach dynamical equilibrium, we propagate the 
noisy wavefunction in time, using a plane- wave basis with 
wavevectors up to k s . Hence, we create macroscopic ther- 
mal fluctuations in the low energy modes, with an ap- 
proximate Bose-Einstein distribution 0, 03. Finally, 
when the mode populations have stabilized, we increase 
the cutoff wavevector to k m , add quantum fluctuations 
to all modes, and begin the merging process 15]. This is 
a finite temperature classical field method 
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Here, we present simulations without fluctuations, in- 
cluding quantum fluctuations only, and including both 
quantum and thermal fluctuations. We compare our re- 
sults for each method to isolate the effect of the mean 
field, the quantum noise, and the thermal noise, on the 
performance of the interferometer. 

We begin by presenting simulations without fluctua- 
tions, taking r = 5 ms (Fig. [TJ. For clarity, Figs. HJb-g) 
only show the central part of the atom cloud, within the 
region enclosed by the dashed rectangle in Fig. Ufa), at 
various times during the merging process. When A = 0, 



the two clouds combine gradually [Figs. QJb,c)] to form 
a final BEC [Fig. [TJd)] with a smooth density profile, 
which is very similar to the ground state of the single- 
well potential. However, the dynamics are much more 
complicated if A = ir. In this case, the phase shift cre- 
ates a density node between the two approaching BECs, 
which evolves into a soliton shown by the central 
white band in Fig. [He). Since this excitation is unsta- 
ble, it bends [Fig. Q|f)] and then decays into vortices 
[Fig.QJg)], via the snake instability 1|| lit ]. 

The addition of a soliton increases the kinetic energy 
of the wavefunction by ~ Hu) x = 7 x 10~ 31 J per parti- 
cle. When the soliton decays, this energy is transformed 
into the kinetic energy of the vortices. Since the BEC 
is in the Thomas-Fermi regime, in the absence of soli- 
tons or vortices the energy available to expand the cloud 
is simply the mean-field energy, which is approximately 
0.5nx4iTh 2 a/m — 15xl0~ 31 J per particle, where a = 2.9 
nm is the s-wave scattering length and m is the mass of a 
23 Na atom. The soliton increases this energy by ~ 50%, 
causing an enhanced spreading of the cloud when the 
BEC is released from the trap, as seen in experiment [l[. 

In the experiments [l[, changing A was found to have 
little effect on the expanded cloud width for either very 
small (< 5 ms) or very large (> 50 ms) values of r. Our 
simulations reveal similar behavior and enable us to iden- 
tify the underlying physical reasons, which are different 
in the low and high r regimes. When t is very small, the 
merging process is dominated by the impact of the colli- 
sion between the BECs, irrespective of the value of A. To 
illustrate this, Fig.[2]shows the merged atom density pro- 
files for t = t = 1 ms and A = (a) and 7r (b), obtained 
from simulations without fluctuations. Both density pro- 
files have modulations in the ir-direction, caused by inter- 
ference between the colliding BECs. In both cases, the 
subsequent expansion of the merged BEC is largely de- 
termined by the momentum of the colliding atoms, which 
masks any enhanced spreading of the cloud resulting from 
the soliton in Fig. [2fb). When t is very large, the BECs 
recombine at low speed. As a result, the merged density 
profile for A = 0, shown in Fig. [2fc) for t = r = 50 ms, 
is similar to the ground state of the single-well potential. 
In this large r regime, the soliton formed when A — tt de- 
cays before the BECs are fully merged, and consequently 
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FIG. 3: Fourier power spectra, F (k z ), obtained from sim- 
ulations without fluctuations at t — r — 5 ms, for A = 
[solid curve in (a)], 0.87T [filled circles in (a)], and n [solid 
curve in (b)]. The vertical dotted lines indicate k z = ±fc c . (c) 
N c /N versus A curves (t = r = 5 ms) obtained from simula- 
tions without fluctuations (solid curve with crosses), including 
quantum fluctuations (dotted curve with circles), and includ- 
ing both quantum and thermal fluctuations at T = T\ (dashed 
curve with downward pointing triangles) and T2 (dot-dashed 
curve with upward pointing triangles). 



few vortices are produced. For example, the density pro- 
file for t = t = 50 ms and A = tt [Fig. [2fd)] is almost 
identical to the equivalent image for A = [Fig. [ZJjc)], 
except that it contains a single vortex. Since this vor- 
tex has a negligible effect on the subsequent expansion of 
the atom cloud, in experiment the width of the expanded 
cloud is independent of A. 

Since the many vortices formed for intermediate t (~ 5 
ms) are described by high momentum plane wave modes, 
we can quantify their production by taking Fourier trans- 
forms of the merged BEC wavefunction. For simplic- 
ity, we integrate the atom density along both the x— 
and y— directions and then calculate the one-dimensional 
Fourier power spectrum, F (k z ) [2fJ. For t = r = 5 ms 
and A = 0, no vortices are present in the BEC [Fig.QJd)] 
and, consequently, F (k z ) [solid curve in Fig. [3(a)] is 
sharply peaked around k z = 0. In contrast, for A = tt, 
F (k z ) is far broader [solid curve in Fig. [3(b)], because 
the BEC contains many vortices [Fig. [1(g)] . 

Our simulations without fluctuations predict sharply 
resonant vortex production when A = n. This can 
be seen from Fig. 03a), which shows that F (k z ) for 
A = 0.87T (filled circles) is almost identical to the cor- 
responding curve for A = (solid curve). To quantify 
this resonant behavior, we consider the number of atoms, 
N c = J ", F [k z ) dk z , with small wavevectors of magni- 
tude \k z \ < k c = 7 x I0 4 m" 1 in the merged BEC [the 
dotted lines in Figs. [3(a) and (b) indicate k z = ±k c ]. 
The quantity N c is close to N for vortex-free BECs with 
smooth density profiles, whose F (k z ) curve is sharply 
peaked within — k c < k z < k c [Fig. [3(a)]. However, 
N c drops rapidly as topological excitations appear and 
broaden F [k z ) [Fig. Eft))]. In Fig. [3(c), we plot N c /N 
as a function of A for merged BECs with t = r = 5 ms. 
For the simulations without fluctuations (solid curve with 



FIG. 4: Atom density profiles in the y = plane (axes inset) 
showing interference between the two BECs after 8 ms time- 
of-flight with A = at T = Ti (a) and T 2 (b). Vertical bar 
shows scale. 



crosses), N c /N decreases rapidly for A > 0.97T. When 
A = 7r, N c /N « I because almost all of the atoms 
have large wavevectors, and therefore travel quickly from 
the trap center when released from the confining poten- 
tial. Consequently, the resonant production of vortices 
should reveal itself as a sudden increase in the width of 
the expanded cloud as A approaches tt. Simulations in- 
cluding quantum fluctuations [dotted curve with circles 
in Fig. [3(c)] show similar resonant behaviour. However, 
the experiments found a far smoother variation of the 
expanded cloud width with A [l[ . We therefore conclude 
that finite temperature phase fluctuations must also be 
considered to explain the experiments. 

To elucidate the effects of finite temperature, we first 
study phase fluctuations by abruptly releasing two BECs 
from the double-well potential at t = and allowing 
them to expand into each other for 8 ms, as in exper- 
iment Q. Phase fluctuations can then be observed as 
kinks in the resulting interference fringes. Figure [4(a) 
shows an interference pattern in the y = plane pro- 
duced by two BECs which, on average, have two atoms 
of extra noise in each low energy mode with k < Jc= . We 
refer to the corresponding temperature T as T\ [15| . Al- 
though the fringes are very clear, their x-position varies 
along the z— direction due to phase fluctuations, which, 
in effect, make A weakly z— dependent. As we raise T 
by adding more noise, the magnitude of the phase fluc- 
tuations increases and their mean wavelength decreases. 
Consequently, the interference pattern obtained for BECs 
with four atoms of extra noise per mode (corresponding 
to T = T2 > Ti) has larger kinks in the fringes, which 
occur over shorter distances along z [Fig. IHJb)]. 

We now calculate N c /N as a function of A for merging 
BECs at T = Ti and T^. Since each individual simula- 
tion is noisy, we repeat the merger five times for each A 
value and average the results. In contrast to our T = 
calculations, our simulations at T = T± reveal the pro- 
duction of vortices for all A > 0.57T. As a result, the 
corresponding N c /N versus A curve [dashed with down- 
ward pointing triangles in Fig. [3(c)] is almost constant 
for A < 0.5-7T, but decreases approximately linearly with 
increasing A > 0.57T. At the higher T = T2, N c /N varies 
more smoothly with A [dot-dashed curve with upward 
pointing triangles in Fig. [3(c)], as observed in experi- 
ment In the limit of very high T, the two BECs break 
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dence of vortex production on A. However, a very precise 
detection of A = tt may be required, as, for example, in 
the "motion detector" shown in Fig. 03a). In this case, 
either the BECs should be very cold, or approximately 
spherical, in order to suppress phase fluctuations. 

We thank W. Ketterle for helpful discussions. This 
work is funded by EPSRC. 



FIG. 5: (a) Atom density profiles in the y = plane (axes in- 
set) of merged BECs [t = t = 5 ms) calculated, without fluc- 
tuations, for A = 4 x 10 4 z (a) and A = 0.057T sin(10 4 7r,z) + tt 
(b). Large dashed rectangles contain enlargements of density 
profiles within smaller rectangles. Horizontal bar shows scale. 



up into multiple quasi-condensates, and their initial rel- 
ative phase at any particular z co-ordinate is essentially 
random. In these circumstances, the number of vortices 
produced is independent of A, and so the interferometer 
no longer functions. 

We can investigate further the effect of a spatially- 
varying phase on the merger by imposing a particular 
z— dependence on A at t = 0. Firstly, we consider a 
linear phase gradient, A = 4 x 10 4 ,z, which means that, 
initially, the two BECs are moving with a relative speed 
v = 4 X lO^h/m — 0.1 mm s _1 in the z— direction. Since 
this speed is so low, the phase gradient has a negligible 
effect on the center-of-mass position of the approaching 
BECs when r = 5 ms. However, the condition for res- 
onant soliton production, A = tt, is only satisfied for 
2 = ±79 fim. As a result, the merged BEC contains vor- 
tices at these positions only [Fig. [5ja)]. Consequently, 
the initial relative velocity v = h/mSz can be deter- 
mined by measuring the distance Sz between vortices in 
the merged cloud. This distance is ~ 300 times larger 
than the change in the atoms' center-of-mass position 
that the initial phase gradient produces over time r, thus 
magnifying the BECs' relative motion. Furthermore, the 
direction of vortex rotation indicates the direction of rel- 
ative motion: positive dA/dz creates anticlockwise rota- 
tion in the y = plane, whilst negative dA/dz creates 
clockwise rotation in the y = plane. Non-linear phase 
gradients can also be studied using this technique. Fig- 
ure E3h) shows a merged cloud at t = r = 5 ms, for a 
sinusoidal phase imprint A = 0.057rsin(10 4 7rz) + tt, like 



that created in a diffraction experiment 21| . Vortices are 
now observed at z = and ±100 /im. 

In summary, we have shown that the decay of a dark 
soliton into vortices is the physical mechanism that en- 
ables the extraction of phase information in recent ex- 
periments on atom chip interferometry [l[ . The design of 
the interferometer can be tailored for different purposes. 
If all values of A are to be measured, the atom clouds 
should be at finite temperature and elongated to enhance 
phase fluctuations, and hence cause a smoother depen- 
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